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IX’l'ROmiCTION 


A simi'.lalion  model  that  roprosoii  i s.  a su  hma  r i no  ’ lor  a hi;4li 

value  laift'c'l  within  a .'-■pcci  I it'd  ope  rat  in;;'  area  it  deserllted  in  inis 
rejK)i'l  . This  model  was  developed  as  an  adjunct  to  the  lo  rriu  1 a l i on  and 

•4 

1 mp  lomenl  a t ion  ol  a computationally  more  el'l'icient  analytical  model. 

The  s i.mu  1 a t ion  model  served  twi,)  purposes.  First,  the  simulation  model 
provided  a validation  of  the  s 1 ;t  I i .s  t i e;il  in|iii1s  used  lor  the  analyi  ieal 
model.  Speci  I'ical  1 y , tlie  siiitulalion  studies  vtilidated  the  apip  1 i t.'ab  i 1 i 1 y 
of  the  analytical  motlel  for  tie  I ermininp-  rtite  of  encounter  between  Mth- 
marine  and  'tirpets.  Seetnid,  results  obitiined  t lirouph  exert:ise  of  the 
simultitittn  model  provided  a eouveiiient  cheek  of  the  rt.’tisttnabl  eness  o I' 
analytical  model  resttlls. 


.1.  M.  Mourtt;  ".Semi -Ma  rkov  Models  of  Seareh  in  the  I’reseiua?  of  Dt't.oys,' 
NWltC  RM-thl,  SKI  Project  1 1)  1 (i- 111 , Contriiet  NOOt)  1 1-7  1 --(t-O  1 1 Stanloitl 
Hesettreh  Tiisl  ilule,  Menlo  Park,  Ca  1 t I o rn  i a;  .fitly  1!)7  I (PNP  l.'l  t IS  1 |,'IFD) 
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GENERAI,  SIMUI^TION  MODEL  CHAUACTEKISTICS 


Thu  ovcraLl  simulation  mod<'l  is  (.'omtiri  so(i  o|'  tlircu  suparal  u tomput  ( i' 
submoilo  1 s . Each  of  tliosc  submodels  is  ;lc.sipno;l  to  ropi'csuni  a spociiic 
sot  of  operational  conditions  which  may  characl  c ri /.c  1 ho  various  pha.scs 
of  submarino  search  for  a In  ifh  value  lar^i'l  (HVT)  in  a specif  ic-d  ai'oa  . 

The  operational  conditions  l•<■■pl■l'scnl  ed  in  these  su'biiiode  i s ineiudi,’  contin- 
uous submarine  .scareh--\vi  th  and  without  false  tai'Ki'ts  (decoys)  itresi’iit  , 
and  subtiiarine  search  empioyini;  the  sprint -drif  t ttietie  withont  ftilse 
tarKots  present  . 


Each  of  these  subtiuulels  is  prop;rammi.d  iit  the  DASTC  algebraic  l;ingnaii;c 
for  processinij;  tiltd  oiterafion  by  a tiliie-sliare  opertil  i iip;  system  i neoeporat  i tip, 
a CDC-ti'iOO  computer,  liecaustt  of  Lime-sliare  system  limitations,  pa  rt  i ■•u  1 ar  1 
with  rt^g'arci  to  output  of  i-wsulls,  those  submodels  are  cons  i rue  t ed  ;-o  ;c  to 
t act i li  t a t e imkdH'iKkai  1 ope  ra  t .i  otis  . 


Those  motlel  t:har  ae  ier  I s t.  ie  s which  c.j'e  common  to  all  (hive  submodels 
are  ditscrlbcd  and  discussed  in  thi.s  secLion;  wiu'i'c'as  ihh'U  1 tir  i I i e s oi  each 
sul)mo(k!l  are  ulscuhsed  in  the  sections  whirh  follow. 

Tin;  tirimtn’y  oitlpnl  of  rstch  submodel  is  the  olttpsed  I iiin  1o  1I\'T 
del.ee  Lion  by  llie  submttrine  senreher  and  variou.^t  ,st  at  1st  ic;t  I tiuaiil  i I i e;-; 
wiiich  char.'.icteri /e  the  distribution  of  the  time 


to  IIVT  (lelet'tion.  Specific 


aspi'i'ls  ('!'  thr  '.'arioui-,  cuimiiun  model  eluu’Lie  Li-r  i s t i es  are  d i seu.ssed  in  t he 
l>a iM|.; raph.s  wiiieii  Col  low. 

11.1  Ope ra  i i ona  I K rame'vot  k 

In  e\ery  ea.>^i',  the  areti  i -.roii  side  rod  to  bi'  eireular  In  shaiJe  with 
a known  t'a^liUH,  U.  While  tin-  radius  of  the  area  may  be  inpuftial  a.s  a 
p a r anil- 1 er  , it  lias  usually  bei-n  a.s.sij^ned  a value  nl  1100  niii  i , 


For  referenei’  jiurpose.s.  the  I’enti-r  of  llu‘  ai'ea  i.s  matle  to  coincide 
with  I lie  oriiiin  of  a carte.sian  coordinate  sy.slem,  orienfeti  so  as  to 
position  the  .\--;ixi.s  Ln  a horizontal  direction.  Initial  and  .siib.iei|uenl 
par  t ic  i|);uil  po.sitions  arc  expressed  in  lerm.s  of  the  underlying  coorilinafe 
sy.s  t cm  . 


The  ihilitil  posUions  of  the  IIVT  and  decoys  (it  includiai)  are 
randomly  ilistributed  within  the  area.  Tliis  i.s  accomplished  through 
the  u.se  of  random  luanbers  uni  to  rut. '.y  distributed  on  the  interval  [u,l] 
t he  I'o  1 low  i niiuine  ” . 

bet  u ;ind  u be  i wo  nuiiilH't's  drawn  random  tv  t rum  a population  ol 
number, s win  eh  i.s  uniformly  d 1 s L ributi.'d  on  I he  interval  i0,lj.  Then, 
example,  the  initial  po.sition  of  the  IIV'T,  1’  , is  determined  a.s; 


Initial 


CO 

P (Ru  , lirru  ) [ Rii  cos  (Utiu  ) , Ru 

eo  1'  1 ;t  1 

posit  ions  for  flic  dt;eoys  aro  delermined  in  a 


sin  (Rmi^)] 
similtir  manne>' 


i n 


f O I' 


utiti/.Inp'i  ol'  course 


different  random  luuilbtu's. 


In  Ui<‘  ca.sc'  of  Uio  HUbmarinc  .''.c aivliui' 


the  mode'  user  may 


an  Option  to  select  either  an  initial  position  inside  the  area  or  a 

position  on  the  boundary  oj'  ilie  area.  In  ti^is  latter  case,  the  initial 

submarine  position,  P , would  be  (u  is  again  a random  number); 

so  3 

I’  = (R,  2^u  ) = cos  (2nu  ),  R sin  (2tiu  )] 

so  3 ■ 3 3 

In  the  fonner  case,  the  procedure  would  be  similar  to  that  described 

above  tor  dete’rmining  the  initial  position  ol  the  HVT . 

For  the  majority  of  the  cases  investigated,  an  initial  submarine 
position  on  the  boundary  of  the  search  area  was  selected.  This  condition 
was  considered  to  more  accurately  describe  the  actual  operational  situations 
which  might  be  oncoc.iitercd . 

2 . 2 Model  Time 

Each  of  the  submodels  employs  a fixed  time-step  form  of  operation  lo 
record  elapsed  program  time.  Thus,  participant  motion  and  the  occurrence, 
non-occurrencc  of  various  events  arc;  evaluated  at  constant,  discrete  points 
in  time.  This  form  of  model  construction  necessitates  a tradeoff  between 
the  size  of  the  time-step  and  the:  concurrent  actual  running  time  of  the 
program.  Currently,  a time-step  of  O.-h  lioi.n’s  i.s  used.  This  interval 
provides  a .sufficiently  tietailed  dcscj’lption  of  the  interactions  of  the 
various  forces,  while  not  requiring  c-xorbilani  progrr.m  running  time. 

Of  course,  the  size  of  the  time-step  m.ay  be  inputted  prior  to  progi-am 
uti  lizal inn . 
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2 . 3 Subin afino/ir\'T  'Decoy  Moti on 


Onco  the  initini  positions  of  tlio  various  units  arc  determined,  the 
;iu5tion  ol  each  is  the  result  of  a pseudo-random  process.  That  is,  tlic 
elioice  of  direction  is  purely  random,  but  the  unit  speed  and  longtli  of  time 
on  each  leg  are  prodcterinined . Further,  because  of  the  finite  nature  of 
th'.’  boundary  of  the  area  as  well  as  other  operational  cons  idorat  i on.s , 
the  motion  of  the  various  units  may  bo  inhibited  to  conform  to  those  con- 
ditions. Area  boundary  effects  are  discussed  in  this  section;  other 
operational  factors  which  inhibit  unit  motion  are  discussed  in  the 
appropriate  sections  which  follow. 

A typical  search  initiation  situation  is  illustrated  in  Figure  2-1. 
Where,  as  previously  described. 


P 

CO 

and  P = 

so 

The  Inilial  direction  of  HVT 
that  is,  the  initial  HVT  movement 
area.  This  move” '^nt  is  indicated 
P in  Figure  2-1.  The  searcher, 


(Ru  , 2nii  ) 

2 1 

(R,  2tiu^) 

motion  is  tlien  the  angle  6 - 2mi  , 

CO  4 

is  adially  away  from  the  center  of 
by  the  dashed  line  originating  from 
tjn  the  other  hand,  moves  initially 


CO 

toward  the  cc”'tor  of  the  area  in  a direction  equal  to  9 = 2nu  . It 


so  a 

.sliould  bo  noted  that  the  angles  Q and  0 ,'ire  measured  countei  c lock - 


the 


CO  so 

wise  from  the  positive  \-axis  and  do  not,  therefore,  correspond  to  tlie 
usual  East  West,  North-bouth  orientation.  Decoys,  if  not  stationary,  move 
initially  J.n  a manner  similar  to  that  exhibited  by  the  HVT.  No  furtlier 


b 


comment  will  be  made  in  this  section  retjarding  subsequent  decoy  movement. 


for  tliis  aspect  of  tlio  model  formulation  will  bo  discussed  in  tlxe  section 
deal  in*;  witli  the  peculiarities  associated  with  search  in  tlie  <iocoy 
field  envi  I'onment  . 

Control  over  tlie  movement  of  the  IIVT  and  tlic  subniarine  is  exercised 
in  two  ways;  by  specifying,  as  a program  input,  the  length  tif  each  track 
segment;  and  automatically,  by  "reflecting"  either  the  II\n'  or  tlie  subma- 
rine off  the  boundary  of  the  ai-ca. 

As  an  i I lus  t i-ation  of  tliese  two  procedures,  consider  that  the  length 

of  the  irirT  t:rack  segmonts  lias  been  specified  (inputted)  as  a maximum  of 

L umi . Then,  using  the  situation  depicted  in  the  previous  diagram,  il' 
c 

Y + L H,  the  lIVT  will  move  along  the  dashed  line  until  either  a distance 
c 

I,  from  P is  reached  or  the  im"  is  detected  by  tlie  submarine, 
r CO 

For  cither  condition;  HVT  completes  movement  of  h nmi  oi’  HVT  at  the 

c 

boundary  of  the  area,  a new  direction  of  movement  is  determined  by  calcu- 
lating a "turning  angle,'  0' , such  that 

0'  - 2mi  (u  a random  number) 

and  the  new  dii'oction  of  movement  is  then 

0 --  0 I O' 

now  old 

subject  to  the  constraint  that,  if  the  IIVT  is  outside  tin;  area,  O mu.st 

new 

result  in  returning  the  HVT  to  a posi  ion  within  tiio  area. 
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It  should  be  apparent  from  the  foregoing  discussion  that  the  speed 


of  the  irVT  (and  the  submarine)  is  specified  prior  to  beginning  the  model 
run  and  is  constant  throughout  the  run. 

Movement  >£  the  submarine  is  accomplished  in  a manner  exactly 
similar  to  that  described  for  the  IIVT . 

It  should  bo  noted  that  the  effect  of  keeping;  the  submarine  within 
the  area  boundary  reflects  an  implicit  a.ssumption  that  the  submarine 
possesses  perfect  information  concerning  the  size  and  location  of  the 
search  area.  If  this  were  not  the  case,  it  would  bo  necessary  to  overiaj^ 
the  circular  area  with  an  area  representative  of  the  degree  of  submarine 
intelligence  information.  Such  a submarine  area  might  bo  larger,  smaller, 
and/or  offset  from  the  actual  area. 

2 , 4 Detection  of  the  HVT  and/or  Decoy 

In  each  of  the  submodels,  the  detection  capability  of  the  submarine 
against  either  the  HVT  or  the  decoy  is  described  by  a definite  range 
probability  law  or  "cookie  cuttor."  That  is,  for  scHtic  range  the 

probability  of  detection  is 

(l  r < 

p [Dot]  I 

(O  r > R 

v.horc  r is  the  range  between  the  submarine  and  either  the  HVT  or  the 
decoy  at  the  'nd  of  a particular  time-stop.  In  Hiis  regard,  it  is  po.ssihle 
for  the  submarine  to  pass  the  HVT  or  diicoy  within  detection  range  during 
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a time-step  without  a detection  being  recorded  if  the  range  between  the 
final  positions  is  sufficiently  large.  Again,  judicious  selection  of  the 
time-step,  keeping  the  possible  relative  speeds  in  mind,  can  serve  to 
miuimizt  this  possibility. 

In  actual  operation,  a value  of  may  be  specified  for  both  the  HVT 
and  the  decoy  which  reflects  characteristics  of  the  HVT  and  decoy,  the 
various  envirorjnent al  effects,  and  the  other  factors  which  determine 
detectability  of  HVT  and  decoy.  Obviously,  one  method  of  reflecting  the 
efficiency  of  the  decoy  in  accurately  impersonating  the  HVT  is  to  make  the 
detection  range  of  the  decoy  approach  that  of  the  HVT. 

2 . 5 Input  Parameters 

In  a slight  departure  from  the  prcvioius  format  of  presenting  only 
the  comn\on  features  of  the  various  submodels,  all  of  the  necessary  input 
parameters  for  each  of  t;he  submodels  are  presented  in  Table  2-1.  The 
entries  in  Table  2-1  are  so  arr.nnged  as  to  present  tho  c(5mmon  inputs 
required  and  then  the  specific  inputs  required  for  each  submoilol.  Foi' 
c.ase  of  refurence,  the  BASIC  parameter  symbology  for  the  various  input 
quantities  is  indicated  in  Table  2-1  rather  than  the  more  tiaditional 
mathematical  notations . 

Bocausc  of  the  random  a.spects  introduced  Into  the  models,  tlu;  ri,;.sults 
from  each  of  tho  submodels  are  amenable  only  to  statistical  intei'prutat i on 
For  this  reason,  it  is  necessary  to  replicate  the  model  result.s  many  times 
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Tubir  a-1 


Number  of  PioK»‘um  lUrplloutlons  (Ni>U‘t  NUbS’LI  ai  OUO  for 
Conti nuuii«  Hearcfi--Wi tl)  rnl.Mo  Submodel) 


inuoUH  Som*cii — Without  Km Irto  TarKct.s 


lUuiiun  of  o.V  UVT  by  SulHmu-liU' 

Submurlne  iJeurelt  Velocity 

Lonjith  of  Nubmui’lne  Trm*h  Suginoiil  CCOUHS  i U) 


ni./tn  1 i(.  Seoreh — Wfthoul  rulsc  TnrK^t^ 


yuhii\taim;  iSprlnt  Velocity  <VSJI»‘Sl’i’  S U) 

Subiimrliio  Drift  VolocJly  (VSDf*l)l'i»  i H) 

Suluiiui'iiuj  Sprint  i’oi'iotl 
S'lbiiiuriuo  Drift  Pi'i'ioU 

lliulius  of  Dotectioii  ot  MVT  by  Suhmnrlne  Durini;  Sprint 
U'-ulAMH  of  Detection  of  HVT  by  Submui'lno  DiirluK  Di'i  f t 

Proj^i'niii  UoHtart  Puiumetera — Uh«?  In  t'n.se  t»f 
Pi-o>truiii  i;xeeuti<iM  laterruptiou.  liilllul  Vuluea: 
ij-‘iij':  . i,  uiw;  0.  XV  . i 


Con  t iiiuoUH  Son  i*rb--Wi  I b P»  I sc  Tu  rt;e  I I'i  c 1 <1 


Subrnwrino  Seurcli  Velori ty 

l.i’iibth  of  Sul)inu»’ine  Trueh  ScKiiieiii  ((!OUUS  II) 
LeuKth  of  Decoy  Tr»ch  [l  l,y,3,‘J,0l 

tcouiuKi)  i h] 

iJiM  ,1-  Vi'iiii'ity  [i  Ivih:(j)  > o] 

D<'vi)y  tUofi.sliiealljjn  Time 
IH’Coy  ‘''t\irmjd  Off”  Timo 

K.uUu.s  i)f  Def.ectitm  of  Doct>y  bv  Sulimarlm? 
Uiuliu.H  of  I.>elneiion  of  IIVT  by  •‘juliinni’l  no 


for  oacli  sot  of  operational  conditions  which  arc  lo  lie  invis  I igal  ed  . 


'rhiis  . 


provision  is  included  in  eacli  submodel  for  user  spec  i I rc  a t i on  of  the  iiumhor 
of  replications  of  model  results  desired.  It  siiould  be  noted  liial  , hecause 
of  tile  increased  eomput ational  time  duo  lo  the  added  presence  of  decoy,  t lu> 
numlicr  of  replications  for  the  continuous  scai'ch  witli  the  decoy  suli- 
model  must  be  limited  to  loss  ttian,  or  at  most,  500, 

All  other  entries  in  Table  2-1  are  either  se  1 f -e.xp  1 anatory  or  ari’ 
discussed  olseivhure  in  this  technical  note. 

2 . G Submodel  Outputs 

In  each  case,  the  principal  uuantity  measured  is  the  elapsed  tiiiu’  to 
first  detection  of  the  IIVT  by  I he  submarine.  Since  dolection  is  charac- 
terized by  a definite  rant'e  probability  lav;,  first  detect  ioh  is  equivalent 
to  first  cnc:ounter,  where  encounter  occur.s  whonevijr  the  rann'e  between  the 
.submarine  and  the  HVT  is  less  than  some  predetenninod  value. 

The  output  of  each  submodel  con.sists  tlieu  cif  the  .sruiiple  mean, 
variance,  and  standard  deviation  of  the  t lm<;  to  fir.st  deteriioii.  I'To'  t lie 
two  continuous  search  sulimodels,  these  values  ari>  indicated  by  the  llA.SIC 
sjmibols:  MllANC,  VAHC , aiul  STllEVC , respectively.  Fin-  the  sprint  drill 

submodel,  the  output  foi'mat  has  been  moilifiod  to  tirint  ou'  rlie  full  titles: 
"STANDAUD  UEVI;\TI0N , ” etc  . 
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In  computini';  thu  yample  variancL-,  t liu  con.si  stinil,  i_’s  t imat.  or  formulation 


lor  the  population  varianoo  is  usoct,  that  is: 


M 

- V (T  - T)- 
M ^ i 


i 1 


/ 1 ^ \ 

where  T is  the  sample  mean  I = / 'i'.j  ^ sample  variance. 

The  output  of  the  two  continuous  searcli  submodels  also  includes  b.ih 
a frequency  count  and  a cumulative  i'rcquency  of  occurrencu  of  first  lie- 
tection  as  a function  of  elapsed  time  interviil.  Further,  the  output,  of 
the  continuous  search  with  llu'  decoy  submodel  includes  a listing  ot 
specific  times  of  first  encounter  of  the  HATF  with  the  subniarine  lor  etich 


of  the  replications. 


3. 


CONTl'INUOUS  SEAHCIl  - NO  DKCOV  SUIiNK)DEL 


This  submodol  reprc'sc'iits  the  basic  strucLui'c  io3'  ail  ol  thi'  sul)~ 
models  constructed  to  date.  As  such,  the  internal  structure  oi  lliis 
submodel  is  exactly  similar  to  that  described  in  tlie  prececlini;'  st  ctiun, 

A single  submarine  searcdiei-  seeks  a single  HVT  within  a delineatc'd 
objective’  opi-rating  area.  Tlu'  two  units  movi’  in  a pseudo-random  manner, 
staying  within  the  objecttvi'  ai-ea  at  all  tiini’s.  Each  case  oi-  replicalion 
is  terminated  at  the  instant  oi  initial  submarine-IIVT  eitcounleiv'detect  i on 
The  results  obtained  f'roni  c>xercising  lliis  submoilel  lor  ranges  oi  Llii.’ 
various  input  parameters  provide  the  baseline  data  lor  evalualJ.ug  the 
elfectiveia.'ss  oi  various  tactical  procedures,  ACM  emp  1 os’iiicnt  policies, 
and  combinations  ol  both. 

A typical  example  oi  the  type  oi  motion  whieh  this  subniodL' 1 generate 
records,  and  evaluates  is  shown  in  Kig.  3-1. 

An  example  oi  tlie  results  obtaiiU'cl  irom  exercising  tliis  sul)modet 
is  presented  in  Fig.  3-3. 


submarine 


Figure  3-2  sample  results  from  the  continu* 


SPRim'/DRIFT  SEARCH  - NO  DECOY  SUIiMODEL 


-•I. 

This  submocU'l  difii'i's  irorn  the  baseline  submociol  in  two  principal 
aspects;  the  manner  in  which  the  motion  of  the  submarine  is  rt'^ulateci 
and  the  manner  in  whieli  detection  of  the  HVT  may  occur. 

As  indicatefi  in  the  title  of  this  submodt- 1 , the  submarine  employs 
a tactic  of  first  executing  a high  spet-d  sprint  and  then  a sitiw  spied 
drift  to  attempt  deteetion  of  the  HVT.  Thus,  motion  of  thi'  submarine  is 
characterizi'd  by  tour  input  parameters;  tlie  sprint  spei-d,  the  sprint 
period,  thi-  drift  speed,  and  the  drift  period.  In  this  manner,  tlie  single 
length  of  submarine  track  segment  (COUUS)  specified  in  the  IjuselLiie  sul)- 
modei  is  replaced  by  two  li'ack  lengths  equal  t.o  (VSSP'SSP)  and  (VSDE  • Dl’P) , 
respectlvi'ly.  (Si.e  Table  U-1  for  definition  of  symbois.)  At  the  end  of 
each  of  these  track  Lcngtlis,  new,  random  lieadings  are  determini'd  for  t.iu' 
submai’ine.  This  procedure  is  illustrated  in  tlu'  example  jiloi  of  sulimarine 
and  HVT  movement  presented  in  Figure  d-l. 

Under  actual  conditions,  a submarine  is  effective  !>■  ” aei>ust ic;i  1 .1  y 
blind"  during  tlie  sprint  period.  Tlieri'  may  be  ins  uuicisi , however,  when 
the  speed  seli'cted  for  the  sprint  tactic  is  not  so  high  as  to  compUti-l) 
eliminate  the  detection  capability  of  the  submarine,  especially  11  the 
submarine  and  the  HVT  sliould  pass  close  alioard  one  another.  For  this 
reason,  prnv  i s i cm  is  included  in  Uiis  subinodi’l  for  speeifying  a snlmiai’iue 
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detection  capability  during  tne  sprint  period  although,  in  all  runs  to 
date,  this  value  has  been  set  so  low  as  to  effectively  preclude  such 
detections. 

lour  programs  are  required  to  exercise  this  submodel.  These  programs 
are:  "SPRINT,”  ''llLECR,”  ''UPDATE,"  and  "STATICS."  The  last  of  these, 

"statics,"  performs  the  statistical  analysis  of  the  results  of  the  various 
replications  and  prints  the  summary  data. 

The  program,  "SPRINT,"  comprises  the  main  program  of  this  submodel 
and  models  all  of  the  details  of  the  search.  As  such,  this  program  requires 
specification  of  all  of  the  initial  input  parameter  values  indicated  in 
Table  2-1.  Prior  to  the  execution  of  SPRINT,  an  indepenpJent  subroutine, 
LITBAL,  must  be  called  into  the  processor.  This  is  accomplished  by  the 
system  command,  "Get,  LITBAL." 

The  two  programs,  "fiLECr"  and  "UPDATE,"  are  utility  progtims  that 
manipulate  perraanoni'.  program  data  files.  These  programs  were  found  to  be 
necessary  in  order  to  safeguard  against  loss  of  ac;cumulated  data  in  the 
case  of  mac)iine  failure.’  or  program  interruption  during  lengthy  program 
executions.  The  three  paraincte.'’s , "LFILE,"  "UIEC,"  and  "ip"  are  used  to 
recover  and  restart  program  execution  should  it  be  interrupted.  (LEIUO 
is  the  BASIC  symbol  for  "Last  file  Updated,"  and  LREC  is  the  symbol  for 
"Last  Rcec;rd  Created.") 

The  output  from  the  SPRINT  program  i.s  printed  on  magnetit;  l.ape  aftir 

This  output  is  tlu'n  available  for  eithei- 
IP 


every  20th  replication. 


statistical  analysis,  using  the  STATIC  program,  and/or  permanent  storage 
on  magnetic  tape  using  the  UPDATE  program. 

An  example  of  the  results  obtained  from  exercising  this  submodel 
is  presented  in  Figure  4-2. 
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5. 


CONTINUOUS  SEARCH  - DECOY  FIELD  SUBMODEL 


As  the  title  implies,  this  submodel  iucorporates  the  etleets  on 
searcher  effectiveness  due  to  the  presence  of  decoys  within  a specified 
area.  This  innovation  i-cquix-es  several  deiiartures  from  the  procedure's 
established  for  the  baseline  submodel. 

As  before,  the  submodel  Includes  a single  submarine  searcher  and  a 
single  HVT,  Due  to  limitations  on  computer  core  storage  in  the  present 
time-share  system,  the  number  of  decoys  must  be  restricted  to  five  or 
less.  Increased  capacity  can  be  achieved  through  the  use  of  utility 
programs  such  as  those  described  in  the  preceding  section,  but  this  has 
not  been  done  at  this  time. 

Either  stationary  or  moving  doeoys  can  be  simulated  through  appro- 
priate specification  of  the  decoy  velocity,  VDC(I).  Because  of  some  ol 
the  computations  involving  VDC,  however,  this  value  cannot  be  zero. 
Therefore,  the  slalionary  doetjy  is  repi’oseuted  by  making  VDC  vi.'ry  sinai], 
albeit,  positive. 

The  HVT  and  the  decoys  are  initially  positioned  within  the  area; 
the  submarine  searcher  oh  the  boundary.  All  of  thi'  units  are  pemiitted 
to  move  in  the  pseudo- random  manner  de.scribed  before  subject  to  the 
provision  that  no  decoy  may  approach  1 e UVT  closer  than  (UlR  + H2K)nmi. 

(See  Tabic  2-1  lor  definition  of  symbols.)  This  co  ■ '.trai  nt  insuri's 
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that  tho  influence  fields  (detection  radii)  of  the  HVT  and  the  decoys 
shall  never  owrlap.  If  the  decoys  are  capable  of  movement,  the  HVT  is 
allowed  to  move  in  an  unrestricted  manner  and  the  decoys  are  required  to 
remain  clear  of  the  HVT.  Obviously,  if  the  d ’Coys  are  considered  to  be? 
stationary,  i.c.,  VDC  is  very  small,  this  condition  must  be  reversed. 

If  overlap  should  occur  at  the  time  of  determination  of  initial  positions, 
the  initial  positions  of  the  decoys  arc  adjusted  s<?  as  to  eliminate  the 
overlap. 

If  the  submarine  searcher  comes  within  detection  ranEc  of  a decoy, 
the  submarine  is  considered  to  spend  a period  equivalent  to  TC  hours 
classifying  tlte  decoy.  During  this  time,  the  submarine  is  precluded  from 
making  new  detections  on  either  the  HVT  or  other  decoys.  The  period  TC 
is  also  sometimes  referred  to  as  "decoy  capture  time  or  decoy  hold 
time, " 

During  tho  period  TC,  tin-  decoy  continues  to  move  in  accordance  with 
the  appropriate  pseudo-random  proeoduros.  The  submarine  searcher  is 
assumed  to  move  in  consonance  with  the  decoy.  The  position  of  the  decoy 
at  the  end  of  the  period  TC  is  taken  to  be  the  re-start  position  of  the 
submarine  searcher. 

To  simulate  submarine  memory  of  classified  decoy  locations,  the 
decoys  are  "turned  off"  by  the  model  for  a time  TM  following  the  classi- 
fication period,  TC . This  permits  the  submarine  to  clear  the  immediate 
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area  of  the  decoy.  Obviously,  the  period  TM  should  lx.'  made  at  least 
long  as  RIR/VS  when  specifying  this  input  value. 

If  it  should  occur  that  the  submarine  is  within  detection  range  of 
two  or  more  decoys  at  any  given  instant,  the  model  assumes  that  the  sub- 
marine will  investigate  and  classify  the  nearest  decoy,  ignoring  thi.' 
others.  If  ever  the  Hl'T  is  within  range,  the  model  assumes  that  it  will 
be  detected  and  classified  by  tbe  submarine  without  regard  for  whatever 
decoys  may  also  be  present.  As  in  all  other  cases,  the  instant  of  first 
detection  of  the  HVT  by  the  stibmarine  terminates  the  specific  replication. 

A typical  example  of  the  motion  generated  lor  the  various  units 
within  this  submodel  is  presunted  in  Figure  5-1.  In  this  example,  4 decoy 
are  present  within  the  objective  area  (Dl,  D2,  D3,  and  DA  in  Figure  5-1) 
along  with  the  HVT  and  submarine  searcher. 

An  example  of  the  type  of  results  obtained  from  exercising  this 


submodel  is  shown  in  Figure  5-2. 


EXPECTED  TIME  TO  N\/T-SUB  ENCOUMTER 


FIGURE  b-2  SAMPLE  RESULT'  FROM  THE  CONTINUOUS  SEARCH  — ULCOY 
FIELD  SUBMODEL 
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NWKC  KM-G'l,  CUT  Prujicl  I OHi-S'l  Ti , (^oiilruo:  NOOUlT-7  1 -c-O]  j p;  S'anl'ord 
Ho.-ii'urch  Ins  ti  Lute,  MonLo  Park,  CaUrorula;  Ju  l y 2'.)7  I (UNCIAdSlPIE))) 
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Appoiulix  A 


HiOGUM'I  LISTING  FOR  CONTINUOUS  SFARCll 
NO  DKCOY  SUUMODRL  (MOO  1) 


4 


MODI 


06/ 1 R/71 


00100  PHOORAM  MODI ( INPUT. OUTPUT) 

00110  DIMENSION  TND(  5)  . XD(  5)  . YD(  5)  . ANOn<  5)  . I DN(  5),  COUKDC  5)  , 
00120+  N0PD(5) 

00130  DIMENSION  VD<5) 

001i)0  DIMENSION  TDCVAC  1000) 

00150  DIMENSION  TDCOYC100) 

00160  DATA  COUKD/5..  5.  . 5..  5.  . 5./ 

00  170  DATA  VD/ 1 o , 1 . . 1 . . 1 . . 1 . / 

00180  R*200. 

00190  NINT=500 
00200  ITOT5=0 
00210  IP1=1 

0 0220  IP2=5 

00230  PRINT. tC* 

002A0  PRINT.  *CVA  RANDOM  START  IN  AREA* 

00250  PRINT.  «SUB  RANDOM  START  ON  POUNDARY* 

00260  PRINT. *C* 

00270  CALL  SECONDCX) 

00280  Y»RANE(X) 

00290  XX=0. 

00300  PRINT. +ENTER  VC.  VS.  CO URC.  CO URS.  RO* 

0 0310  READ.  VC.  VS.COURC.  COURS.  RO 
00320  1 = 0 

00330  PRINT. ♦enter  NREPLI* 

003A0  READ.NREPL! 

00350  TIME»0. 

00360  HRS=0. 

00370  N DECOYS 0 

00380  TDELAYs  (-9.*VS)/50.-*-5.9 

00390  KOUNTS0 

0 0A00  KONDS0 

00A10  KONC=0 

0 0A20  PRINT.  ♦ CO  UN  TER.  T ELASPED.T  ENCO  UNTEKF' D* 

00A30  202  CONTINUE 

004A0  INC=0 
0 0 450  INS=0 
0OA60  TNC=0. 

00470  TNS=0. 

00480  HR5=0. 

00490  nFA=S7.+( 17./60- ) 

00500  TD=0. 

r*m(*  «r7i  rn 

Yi  rj  ^ I V)  iv U r u “♦  frj 

00520  NOPS=0 

00530  DO  2 I=1,NDFC0Y 

00540  NOPD(I)=0 

0 0 550  1 DNCI)  = 0 

00560  TND(I)=0, 

00570  2 CONTINUE 

00580  200  CONTINUE 

00590  ANGC=6.28319*RANF<XX) 


A-2 


MODI 


06/10/71 


00600 
00610 
00620 
0 0630 
0 0660 
0 06S0 
00660 
00670 
0 0680 
0 0 690 
0 0700 
00710 
00720 
0 0730 
0 0 740 
0 0 7S0 
0 0760 
00770 
0 O780-I- 
00790 
0 0800 
0 0810 
0 0820  + 
00830 
0 0840 
0 08  50 
00660 
0 0870 
0 0880 
00890  + 
0 0900 
0 09  1 0 
0 0920 
00930 
00940 
0 0950 
00960 
00970 
0 0980 
00990 
0 1000 
0 1010 
0 1020 
0 10  30 
0 1040 

0 1 0S0 
0 1060 
0 1070 
0 i 080 
0 109  0 


RAnUS=R*RANFCXX) 

XD=  RADI  1 5*  CO -SCAN  GO 
YC=RADUS*SIN(  ANGO 
ANGS=  A. 28319*  RAN F<  XX) 

RAD=K 

XS--RAD*C0S(  ANGS) 

YS=HAP*  ?1N( ANGS) 

I FCNDFCOY. FQ. 0)  GO  TO  800 
DO  4 I=l,NDFCOY 

CALL  INDC0Y(  ANGD(  I ) , X D(  I ) , YD(  I)  ^ X C#  YC>  X S,  Y 5^  XX,  R,  VC*  > 
TNDC I )=  TNDC I ) * I - 
4 CONTINUE 
800  CONTINUF 
DT=0. 5 

15  CONTINUE 

TIMEaTIME+DT 

HR5=HRS+DT 

CALL  LI  TBAL(  TWC*  XC*  YC,  ANGC*  XNOW*  YNOW.  INC*  VC*  COLIKC.  XX* 
NOHC*  TDFLAY*  TO 
XC=XN0W 
YC=YN0 W 

CALL  LITBAL(TNS*X5*  YS*ANG5*XN0V-*YN0W,  IN5*  VS.COUKS.XX* 
NOPS*  TDFLAY*  TO) 

XS=XN0V,' 

Y5=YN0V 

I FCNOHS.  F(J.  2)  GO  TO  360 
I FCNDECOY. FQ. 0)  GO  TO  350 
DO  350  l=l*NDECOY 

CALI.  LI  THAL(  TND(  I ) * Xn<  1 ) . YD(  I )*  ANGDC  1 ) * XNOW*  YNO  W,  I UNCI  )* 
vnc  I ) , cmiKDC  i ), xx* nopdc  i ) * tdflay,  tu) 
xn<  I ) = XN0W' 

YDC  I )i:YNOW' 

3 50  CONTINUF 

CALL  MFFTC  VC,  VS,  X C*  Y C*  X S,  Y S*  I EN  ) 

IF<IFN.FQ«1)  GO  TO  100 
GO  TO  300 
100  CONTINUF 
KOUNT=KOUNT+  1 
KONC=KONO  1 
TUCVA<KONC)sHRS 
I PRIN T= I PI  * 1 0 

IFC  I PRINT.NE.KOUNT)  G^  to  1000 
I P 1 ~ I P 1 + 1 

PRINT  95*K0UNT,HK5*  time 
95  FORMATC I 1 0, 2F 1 1 . 4/ ) 

1000  CONTINUF 

I FCKOUNT.  F'Q.NKFPLl  ) GO  TO  9 7 
1 3TAT= I P2*  1 0 

I F( I STAT.NF.KOUNT)  go  to  1010 

CAI.L  T!  C(  TOCVA*  TDCOY*KONC*KONU.  EXFl*  VAR*  I SIAT,  VS) 
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R1100  C LL  FHE0(  TDCVA,  KONC*  I STAT.  NRFPLI  * NINT,  I TnT5) 

0 1110  PRINT,  *TIMF  RT  5000* 

0 1 120  PRINT,  I TRT5 

0 11 30  I P2= I PS+5 

0 n 40  GO  TO  202 

0 1 1 50  1010  CONTINUF 

0 11 A0  GO  TO  202 

0 1170  300  CONTINUE 

01180  I FCNOECOY.EG.O)  GO  TO  15 

01190  no  361  I=1,NDFC0Y 

01200  CALL  NEET(  VC,  VS,  XD<  I ),  YD<  I ) ,XS,  YS,  I EN> 

0 12  10  IFCIFN.FQ.n  GO  TO  500 
01220  361  CONTINUE 

0 1230  I FCKOUNT.  EW.NKEPLn  GO  TO  97 
0 1240  GO  TO  15 
01250  310  CONTINUE 

0 1260  IF(NLPS.FQ.I)  go  TO  15 
01270  PRINT, *MFET  DECOY* 

012H0  PRINT;  I,XD(I  ),  YD<  n 
0 1290  PRINT,XS,  YS 
01300  1KFFP=I 
01310  NOPS= 1 
0 1320  N0PDCn=2 
01330  TD-0. 

0 13  40  GO  TO  1 5 
01350  360  CONTINUE 

01360  NOPS=0 
01370  NOPDC IKFEP)=0 
01380  XS=XncIKEEP) 

01390  YSaYnCIKFFP) 

01400  DO  370  I=1,NDEC0Y 

01410  I F( I . FO. IKFFPl  GO  TO  370 

01420  CALL  MFFT<  VC,  VS,  XD(  I),  YDU  ),  XS,  YS,  I FN) 

0 1430  IFCIFN.FQ.n  GO  TO  310 
01440  370  CONTINUE 

01450  IKFFP=0 
0 1 460  GO  TO  15 
01470  99  CONTINUE 

0 1480  500  PRINT  ,*DFCOY  ENCOUNTERED,  T ENCOUNTERED* 

0 1490  PRINT,  i.HRS 

01500  98  CONTINUE 

01510  KOND=KOND+ 1 

01520  TDCOY(KOND)=HRS 

0 1530  KOUNT-KOUNT*  1 

01540  GO  TO  202 

01550  97  CONTINUE 

01560  I FCKOND. EQ. 0)  GO  TO  504 

0 1570  504  CONTINUE 

01580  CALL  TIC(TDCVA,TnCOY,KONC,KOND,FXPT, VAR,NREPL1, VS) 
0 1590  CALL  FRFQ(T[)CVA,KONC,  ISTAT,NREPLI,MINT,  1 TGT53 
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0 1fi00  PRINTS  *FRFQUE:NCL'  GT  5000* 
01610  PRINT,ITGT5 
0 1 620  END 


M ODl 
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0 1630  SLIRHOUTINF:  STATC  TIMF.  I gout,  FEAN,  STDFV,  VAR) 

0 16-50  DIMENSION  TIME(1000) 

01650  SUM=0. 

01660  DO  5 I=l,ICOUT 
0 1670  SIIM=SI)M+T1ME(I  ) 

016R0  5 CONTINUE 

0 1690  FEAN=  SL'M/FLOAT<  I GOUT) 

0 1 700  VAR=0. 

01710  DO  10  I=l.ICOUT 

01720  VAR=VARt-C  TIMF<  I )-FFAN)**2 

0 1730  10  CONTINUE  • 

0 17A0  VAR=  VAR/FLOATU  GOUT) 

01750  STDFV= SQKT( VAR) 

01760  RETURN 
0 1 770  END 
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0 1780  SURROU  INF  T1C(  TDCUA,  TDCOY,KONC.  .^OND^  FXPT,  VAR,  NHFPLI  , VS) 
01790  DIMENSION  TDCVAC 1 000 ) , TDCOY C 1 00 ) 

01B00  PRINT, *NO  OF  RFPL I CATI ONS+ 

0 1810  PRINT,  NREPL! 

0»820  TDLAY=(“9.'t'VS)/50.  + 5.9 

01830  I FtKOND.  EO.  0)  GO  TO  10 

01840  CALL  STAT(TDCOY,KOND,MEAND,STDFVD,VARD) 

0 1850  PRINT,  MFAND,  STDFVD,  VARD 
0 1860  10  CONTINUE 

0 1870  PRINT, ♦MEANC, STDEVC,  VARC* 

01880  CALL  STAT< TDCVA,KONC, FEANC, STDFVC, VARC) 

0 1890  PRINT, FEANC, STDFVC,  VARC 
01900  I F<KOND. EQ. 0)  RETURN 

01910  EXPT=(P/C l.-P))*(TDLAY+MEAND>+TDLAY+MFANC 
01920  TERM l=P/( C 1 .-P)**g) 

01930  T~RM2= < TDLAY+MEAND) 

0 1940  TERM3=CP/C  1 .-P)  lif^VARD 
01950  PRINT, TERMl , TERM2, TERM3 
0 1960  VAR= TERM  1 »TFRM2+TFRM3+ VARC 
0 1970  RETURN 
0 1980  END 
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0 1990  SUBROUTINF  LITBAL(  TN*  XLASTj  YLAST>  ALAST»  XNO  VNO  IND>  M. 
0 2000+  COURSFj.XX,NOP^  TDFLAY,  TD) 


0 2010 

1F<N0P.EQ.  1 ) 

GO 

TO 

30 

0 2020 

DT=  . 5 

02030 

R-200. 

02040 

I F< IND. EQ.0> 

GO 

TO 

20 

0 2050  FLEN- TN4<DT*  V 

02060  iFfFLEN. GT. COURSE)  GO  TO  20 
02070  DX=(DT*V)*COS<ALAST) 

02080  DY=(DT*V)*S1N<ALAST) 

0 209  0 RNOW=SQRT( < DX+XLAST3 * * 2+ C DY  + YL AST ) ♦* 2) 

02100  I F(  RNOW.  GT.  R)  GO  TO  20 

02110  XNOW=XLAST+DX 

0 2120  YNOW=YLAST+DY 

02  130  GO  TO  10 

02140  30  CONTINUE 

02150  XNOV'=XLAST 

02160  YNOW=YLAST 

02  170  I P(  TD.LE.  TDELAY)  GO  TO  15 
02180  PRINTj*LAST  TO* 

02190  PRINTjTD 
0 2200  IEN=0 
0 2910  N0P=2 

02220  TDa0. 

02230  ALAST=6.28319tRANE(XX) 

02240  IND=1 
0 22  50  TN=1.. 

02260  PRINT>*T  DELAY  END* 

0 22  70  GO  TO  15 
022R0  15  TD-TD+DT 

02290  PRINT, *TD* 

02300  PRINT,  TD 
02310  GO  TO  25 
02320  20  CONTINUE 

0 2330  IND=1 
0 2340  TN=  1 . 

02350  ALPHA=6.28319*RANF(XX) 

02360  PHY=ALPHA+ALAST 
02370  PHY=AM0D(PHY,6.283I9) 

02380  DX=(DT*V)*COS(PHY) 

02390  DY=(DT*V)*S1N(PHY) 

02400  RNOW=SORT( ( DX+XLAST) **2+ ( DY+ YLAST) **2) 
0 2410  I F(  RNOV.'.  GT.  R3  GO  TO  20 
02420  ALAST=PHY 
0 2430  XNOF=XLAST+DX 
02440  YNOW^  YLAST+DY 
0 9450  10  TN=TN+1. 

02460  25  CONTINUE 

0 2470  KFTURN 
0 2480  FND 
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Pig490  SUBROUTINE  INDCO  Y ( AN  GD  1 . XD  1 . YD  1 » XC.  YCj  X YS,XX»  VC^  VS^ 

02500  5 CONTINUE 

02510  ANGD1= 6. 283 19*RANE( XX ) 

02520  XDl  = R*COS(  ANGDl  > 

02530  YD1  = R'*SIN<  ANGDl  ) 

02540  CALL  MEET(  VC,  VS,XniiYDl.XS.YS.IEN) 

0 2550  IF<IEN.F.G-n  GO  TO  5 
02560  RETURN 
02570  END 


MODI  06/J0/7!. 


02580  SUBftOUTlNE  MEEK  VC#  VS,  XC>  YC,  XS,  YS,  I EN  ) 
02590  R0=60. 

0 2 600  RANGE=SGRT(  <XC"XS)**2+CYC-YS>0’t‘2> 

02610  IF(RANGE.GT.R0)  GO  TO  10 

0 2620  IEN=1 

0 2 630  GO  TO  20 

02640  10  IEN=0 

0 2650  20  CONTINUE 

0 2660  RETURN 

02670  END 
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0 2 680 
0 2690 
0 2700 
02710 
0 2720 
02730 
02740 
02750 
0 2760 
02770 
02780 
02790 
0 2800 
0 2810 
0 2820 
02830 
0 28  40 
02850 
0 2860 
02870 
0 2880 
02890 
0 2900 
0 29  10 
02920 
0 29  30 
0 29  40 
0 29  50 
0 2960 
0 2970 
02980 
0 2990 
0 3000 
0 3010 
0 3020 
0 3030 
0 30  40 
0 3050 
0 3060 
0 3070 
0 3080 
0 3090 
0 3100 


SUBROUTINF  FRFQ( TDCVA, KONG# 1 STAT, NKEPL I * N I N T, I TGT5) 
DIMFNSION  TDCVAC  1000) j I FRF<  500) , TCOMC  500) 

DIMENSION  IFREC(500) 

ITGT5»0 

ifostat.ltdNRepli)  go  to  5 
1 SUM=NREPLI 
GO  TO  6 

5 ISUM=ISTAT 

6 TINT*  10. 

TMAX*0. 

DO  7 - 1 , N I N T 

1FRE< J)=0 
IFREC< J)»0 

7 CONTINUE 

DO  10  I » 1 , I SUM 

DO  20  J=1,N1NT 

TCOMPA=FLOAT( J)*TIN  T 

iFrrocvAf  I ) . gt.  5000. ) go  to  40 

IF< TDCVACI ) . GT. TCOMPA)  GO  TO  20 

IFRE(J)*IFRECJ)-*-l 

IFCTCOMPA.lt. TMAX)  GO  TO  10 

TMAXsTCOMPA 

GO  TO  10 

40  ITGT5-'I  TGT5-*-! 

GO  TO  1 0 
20  CONTINUE 

10  CONTINUE 

1P=  1MAX/TINT+  1 
DO  J0  1=1,  IP 
TCOMC I )=FL0ATC I )*TINT 
30  CONTINUE 

IFRFCC 1 )=IFRFC 1 ) 

DO  110  I --2,  IP 

IFRFCC I )= I FRFC I )+ I FRECC I- 1 ) 

110  CONTINUE 

PRINT, * TIME  INTERVAL, FKEGUENCF* 

DO  100  1=1, IP 

I F<  ! FREC  I ) . EQ.  0)  GO  TO  100 
PRINT  35,  TCOMC  I >,  I FREC  D,  I FRECC  1 ) 

35  FORMATCIH  , 2X, F I 1 . 2, 2I  1 0) 

100  CONTI NUF 

RETURN 

END 
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PROGRAM  LIBTING  FOR  SPRINl’/DRIFT  SEAHCIZ 
HO  DECX)Y  SUBMODEL  (SPRINT) 


D--1 


SPRINT 
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0010M 

00110 

00120 

0 0130 
00140 

001  50 
0 0160 
0 01  70 
0 0180 
00190 
00200 
00210 
00220 
00230 
0 0240 
00250 
0 0260 
00270 
002P0 
00? 

0 0300 
0 0310 
00320 
00330 
00340 
00350 
0 0360 
00370 
00380 
00390 
0 0400 
00410 
0 0 420 
00430 
0 0 440 
00450 
0 0 460 
00470 
0 0480 
0 0490 
0 0 500 
0 0 510 
0 0520 
00530 
0 0540 
0 0550 
0 0560 
0 0570 
00580 
0 0590 


PROGRAM  SPRINTC  I NPUT,  D UTPUT,  TAPF  I , TAPE:?,  TAPF3) 

CALL  RFTR( ! , 5HTAPE1) 

CALL  RETR( 2, 5HTAPE2) 

REWIND  3 
USE(LI TBALl 

PRINT, +LAST  FILE  UPDATE?* 

READ,LF1LF 

PRINT, +LAST  REC  CREATED?* 

READ,LREC 

PRINT,* VSSP,  RDS,  SPP?  + 

READ, VSSP,  RDS,  SPP 
PRINT,*VEDF,RDD, DFP7* 

HEAD, VSDF,  RDD,  DFP 

PRINT,*SUB  RANDOM  START  IN  AREA=l,ON  BOUNDARYs 2,  WH 1 CH ? ♦ 

READ. I START 

R*200. 

PRINT, *NO  OF  REPLICATIONS?* 

READ,NREPLI 
PRINT, *VC,COURC?* 

READ, VC,COURC 
XX=0. 

PRINT,*  ENTER  IP* 

READ, IP 
ICOPYsI P*20 
K = LFIL.F 
KOUNT=LREC 
CALL  SECCND(X) 

Y*RANF<X) 

220  CONTINUE 
ANGC=6.28319*RANE-(XX) 

RADUSsR*RANF(XX) 

XC=RADUS*COS(  ANOC) 

YC  = RADUS*SIN<  ANGC) 

ANGS=6.  28319*RANF(XXI 
GO  TO  ( 10,20), I START 
10  RADIJS=R*RANF(  XX) 

GO  TO  25 
20  RADUS=R 
25  CONTINUE 
XS=RADUS*COS(  ANGS) 

YS=RADUS*SIN( ANGS) 

DT=.  5 
TIME=0. 

INC  = 0 
INS=0 
TNC=0. 

TN  S=  0 . 

Y=RANF(XX) 

ISFSW=0 

IFCY.GT..5)  ISPSW=1 
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00600 
00610+ 
0 0620 
00630 
0 06  40 
00650 
00660 
00670 
00680 
00690 
0 0 700 
00710 
00720 
00730 
0 0740 
0 0750 
00760 
0 0770 
0 0780+ 
0 0790 
0 0800 
0 0810 
0 0820 
00830 
00840 
0 0850 
0 08  60 
0 0870 
00880 
00890 
0 0900 
0 0910 
0 0920 
0 09  30 
0 09  40 
0 09  50 
0 09  60 
00970 
0 0980 
00990 
0 1000 
0 1010 
0 1020 
0 1030  + 
0 1040 
0 1050 
0 1060 
0 1070 

0 1 i'j  i0 

0 1090 


CALL  LITBAL<TNC>XC>YC.ANGC>XNOW,YNOW,  INC*  VC*  COURC*  XX 
DT*  R) 

XC=XN0W 

YC=:YN0W 

IF( I SPSW.BQ. 1 > GO  TO  30 

TMSP=TIMF+DFP 

VS=V5DF 

COURS=VSDF*bFP 

RDsRDD 

GO  TO  3 5 

30  TMSP=T1MF+SPP 

VS=VSSP 

COLIRS=VSSP*SSP 
RD=RnS 
35  INS=0 
TNSa0. 

150  CONTINUE 

CALL  MOVS<  TNS*  XS.  YS.  ANGS,  XNOW.  YNOU.  INS*  VS*  COURS,  XX, 

DT,R1 

XS=XNOW 

YS=YNOW 

DI  ST=SQRT(  <XC-XS)**2+(  YC-YS)i"*2) 

IFCDIST.LF.RD)  GO  TO  100 

IFCTMSP.FO.TIME)  GO  TO  110 
IF(TMSP.GT.TIMF)  GO  TO  120 
PRINT,*ERR1- TMSP.LE. TIME* 

110  CONTINUE 

I E< I SPSW. EG. n GO  TO  130 

I SPSW=  1 

TMSP=T!  1E+SPP 

VS=VSSP 

RD=KDS 

GO  TO  1 40 

130  CONTINUE 

I SPSV,'--0 

TMSP-TIME+t)FP 

VS= VSDF 

COURS-VSDF+OFP 

RDSRDD 

S40  1NS=0 

TNS=W. 

12  time=tihe+ot 

C/  L LITBALI  TN(',XC*YC.  ANGC*XNOW*YNOW,  INC*  VC*  COURC,  XX 
DT*  R) 

XC=XNOV; 

YC=YNO W 
GO  TO  1 50 
100  CONTINUE 
KnUNT=KOUNT+  I 
WRITFC3,I60)  KOUNT*'  ME 
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f»  1 1 !30 
0 1110 
0 1 120 
0 1 1 30 
0 I M0 
0 1 1 S0 
0 1 1 #^0 
0 1 1 70 
0 1180 
0 1190 
0 1 200 
0 1210 
0 1 220 
0 1230 
0 1240 
0 1250 
0 1 2(^0 
0 1 270 
0 1280 
0 1 290 
0 1 300 
0 1310 
0 I 320 
0 1 330 
0 I 340 
0 1 350 
0 1 360 
0 1 370 
0 1 380 
0 1 390 
0 1 400 
0 1410 
0 1 420 
0 1 430 
0 1 440 
0 I 450 
0 1 460 
0 1 470 
0 1 48  0 
0 I 490 
0 1 500 

I C 1 fX 

I»  I wj  1 t.' 

0 1 520 
0 1 530 
0 1 540 
0 1 550 
0 1 560 
0 1 570 
0 1 580 
0 1 590 


160  FORMAK  1 10,  F20.  41 
I FCKOUNT.LT.  I COPY)  GO  TO  220 
I F( I PFC.NF. 0)  GO  TO  404 
RFWirin  1 
RF.WIND  3 

WRITF(1,4131  KOUNT 
407  READ(3,1601  I.TIMF 

I F< I . FG.KOUNT)  GO  TO  406 
WRITF(I,160)  I, TIME 
GO  TO  407 

406  WRITE! 1,160)  I, TIME 
ENOFILF  1 
REWIND  1 

CALL  RFPLC 1, 5HTAPF1 ) 

REWIND  3 
I RFOKOUNT 
K=  1 

GO  TO  310 

404  CONTINUE 

GO  TO  (400,410)  K 

400  CONTINUE 

REWIND  1 

REWIND  2 

REWIND  3 

V'RITE<2,413)  KOUNT 
KEAD(1,413)  NREC 
413  FORMAT! 110) 

I REC=0 

I F(NREC»EO.0)  GO  TO  600 
4 1 5 READ( 1,160)  I , TIM  E 
1 FC I . EO.NHEC)  GO  TO  411 
I RFC= I REC+ 1 

WRITE(2,160)  I RFC, TIME 

GO  TO  415 

411  IRFC=IREC+1 

WRITF(2,I60)  I RFC, TIME 

600  RFADC3, 160)  I, TIME 

IFCIREC.EW.KOUNT)  GO  TO  420 

I F( I . EG. KOUNT)  GO  TO  420 

IREC=IRECfl 

WRITE(2,  160)  I RFC, TIME 

GO  TO  600 

420  IRFC-IREC+1 

WRITE! 2,  160)  1 RFC, TIME 

ENDFILF  2 

REWIND  2 

CALL.  REPL!2,  5HTAPF2) 

K = 2 

REWIND  1 
REWIND  2 
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01600  REWIND  3 
01610  GO  TO  310 
0 1 620  410  CONTINUE 

01630  REWIND  1 
0 1 640  REWIND  2 
0 1 650  REWIND  3 
0 1660  WRITE!  1.413)  KOUNT 
01670  READC2.413)  WREC 

0 1 680  I REC  = 0 

0 1690  I F<NREC.  EQ.  0)  GO  TO  605 

01700  610  RFADCg,160)  I. TIME 

01710  IE< I . EQ.NRFC)  GO  TO  615 

0 1720  5REC--IREC+1 

01730  WRITE! 1.160)  IREC, TIME 

0 1 740  GO  TO  610 

0 1 750  615  I KEC=*I  REC+  1 

01760  WRITE! 1,1 60)  IREC, TIME 

01770  605  READ!3,160)  I. TIME 

01780  IF! I REC. EO.KOUNT)  GO  TO  620 

01790  I E! 1 . EG.KOUNT)  GO  TO  620 

01800  IREC=!REC+1 

0 1810  WRITE! 1.  1 60)  IREC. TIME 

01820  GO  TO  605 

01830  620  IREC=IRFC+1 

01840  WRITE! 1.160)  IREC. TIME 

01850  ENOEILE  1 

0 1860  REWIND  1 

0 1870  CALL  REEL!  1. 5H TAPE  1 ) 

01880  K=1 

01890  REWIND  1 
0 1900  REWIND  2 
0 19  10  REWIND  3 
0 1920  GO  TO  310 

01930  310  PRINT, *LAST  REC  CREATED, LAST  FILE  CHEATED* 

01940  PRINT, KOUNT. K 

0 1950  I EIKOUNT.  EO.NREPLl  ) GO  TO  300 

0 1960  1P=IP-M 

01970  PRINT, +MEXT  IP  IS* 

0 1980  PRINT.  IP 
0 1990  ICOPYt=IP+20 
0 2000  GO  TO  220 
02010  300  CONTINUE 

02020  1P=IP+1 

0 2030  PRINT.  *NFXT  IP  IS* 

02040  PRINT,  IP 
0 20  50  END 
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0 2060 
02070-^ 
0 2080 
02090 
02100 
0 2110 
02120 
0 2 1 30 
02140 
0 2 I 50 
02  1 60 
02170 
02180 
02190 
0 2200 
022  i0 
02220 
02230 
0 2240 
02250 
02260 
02270 
02280 
02290 
02300 
02310 


SUBROUTINE  M0VS(TN»XLAST»YLAST»  ALAST.  XNOW*  YNO  Wj  I ND» 

COURSE;.  XX,  DT,  R) 
lFfIND.EQ.0)  PO  TO  20 
nX=m>V>(<COS(  ALAST) 

DY=DTx>  V* SIN ( ALAST) 

RNOW=  SQRTC  ( DX-*-XLAST ) >*‘’*■2+  < DY  + YLAST)  ♦*2> 

I E<  RNOW.  GT.  R)  GO  TO  20 

XNOW=XLAST+ DX 

YNOW=YLAST+DY 

GO  TO  10 

20  CONTINUE 

IND=  1 

ALPHA=6.2R319*RANE(XX) 

PHY=ALPHA+ALAST 

PHY  = AMOD( PHYj  6. 28319) 

DX=DT*V+COS(PHY) 
nY=nT*V*SIN<  PHY) 

RN0W--SQRT<  (DX  + XLAST)0*2-*-(  0Y4-YLAST)»^2> 

1 F(  RNOW.  GT.  R)  GO  TO  20 

ALAST=PHY 

XNOWaXLAST+DX 

YNOW=YLAST-^DY 

10  CONTINUE 

TNbTN+  1 . 

RETURN 

FNP 
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0 01«0+  BUPROUTINF  LI  TBALC  TN  , XLAST,  YLAST»  ALAST^  XNOW,  YNOW,  IND^  Vi.  COUKiiE,  AA 
00120  IF(IND.EQ.0>  AO  TO  20 
0 01  30  FLFN=TN  + DT*V 

001A0  IFCFLFN.GT. COURSE-)  GO  TO  20 
00150  DX=cnT*V)+COS(ALAST) 

001<E0  DY=(nT*V)>!<SlNf  ALAST) 

0 01  70  RN01<'=  SORTC  < OX  + XLA  ST ) * * 2+  ( DY  + YLAST)  ♦ >•>2  ) 

0 01R0  I Ff  KNOW.  GT.  R)  G50  TO  20 
0 0 190  XNOW=XL,AST-^[)X 
0 0200  YNOWsYLAST+DY 
0 0210  GO  TO  10 
0 0220  20  CONTI  NUF 

0 02  30  IND-1 
0 0240  TN=1, 

00250  ALPHA=6.  28319  + KANF(  XX) 

0 0260  PHY-ALPHA+ALAST 
00270  PHY=AM0D(PHY, 6.2H319) 

0 0280  nX=<m>V)*COS<PHY) 

00290  nY=(nT+V)«SlN<PHY) 

0 0300  RNOW=  SOHT<  < DX  + XLAST ) + * 2+  < DY  + YLAST)*  + 2) 

0 0310  I F<  KNOW.  GT.  R)  GO  TO  20 
00320  ALAS1=PHY 
0 0330  XNOV:=XLAST+nX 
0 0340  YNOP=YLAST+nY 
0 0350  10  TN=TN+1. 

00360  RETURN 
0 03  70  FNO 
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00100  PROGRAM  STATICSC  INPI)T,mJTPt.lT,TAPFA,  TAPES, 

0 0110+  TAPE 6, TAPE 7, TAPES, TAPE9, TAPE 10, TAPEl 1 , TAPE  12) 

00120  13  CONTINUE 

00130  PRINT, ♦ENTER  PILE  N0+ 

00140  KFAD,K 

001S0  GO  TO  C3,  3,  3,  4,  S,  6,  7,8,9,  10,  S 1,  12),K 
0 0160  3 CONTINUE 

00170  PRINT, +ERR- FILE  NOT  USE  .TRY  A6AIN  + 

0 0100  GO  TO  13 

00190  A NAME=5HTAPE4 
0 0200  GO  TO  I A 

00210  5 NAME=5HTAPE5 

0 0220  GO  TO  14 

0 0230  6 NAME-5HTAPE6 

0 02  40  GO  TO  14 

0 02S0  7 NAME>5HTAPE7 

0 02  60  GO  TO  14 

0 0270  8 NAME=SHTAPE8 

0 0280  GO  TO  14 

0 0290  9 NAME=5HTAPE9 

0 0300  GO  TO  14 

00310  10  NAME=6HTAPF10 

0 0 320  GO  TO  14 
00330  n NAME=6HTAPE  1 1 
0 0340  GO  TO  14 
0 0350  12  NAMES  6HTAPE 12 

0 0360  GO  TO  14 
0 0370  14  CONTINUE 

003R0  CALL  RETHCK.NAMr) 

00390  REWIND  K 
00400  READ(K,  55)K0UNT 
0 0410  55  FORMATC  110) 

00420  1 FCKOlINT.  NE.  0)  GO  TO  20 

00430  PRINT, +0  REC  ON  FILE  K,TRY  AGAIN  + 

0 0 440  GO  10  13 

0 0 4S0  20  CONTINUE 

0 0 460  SlJM  = 0. 

00470  50  READ(K,30)  I,  TIME 

00480  30  FORMATC  I 10, E20.  4) 

0 0490  SUM--SUM+TIME 

0 0 500  I FC  I * NE  fKOUNT)  GO  TO  50 

00510  ISUM^I 

0 0520  I FCl  SUM.  EQ.KOLINT)  GO  TO  60 
0 0530  PRINT,  ♦FRROR-KOUNT  NOT  MATCH  LAST  RFC  NO^ 

0 0540  60  PRINT,  ♦ I SUM,  SUM* 

0 0550  PRINT,  I SUM,  SUM 
0 0560  FFANa  SLiM/FLOATC  I SUM) 

0 0570  REWIND  K 
0 0 5H0  SlJM«0,. 

00590  80  READCK,30)  I, TIME 
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0 0 600  X=(  TIME- FEAN)  * + 2 

00610  SliM=Sl)M+X 

0 0620  I F(  I .NF.KOUNT)  GO  TO  8 0 

0 0 630  GO  TO  70 

00640  70  VAK=S11M/E  lATv'ISlJM) 

00650  STAN= SQRTC VAR 
00660  PRINT*  ^SAMPLE'  SIZE  I 
00670  PRINT.  ISUM 

00680  PRINT* +MEAN. STANOARD  DEVI  ATI  ON, VARI ANCE* 
00690  PRINT, FEAN, STAN, VAR 
0 0700  RETURN 
0 0710  END 


UPDATE  06/10/71 


00100 
0 0110  + 
00120 
0 0130 
0 01  40 
0 0150 
0 01  60 
0 01  70 
0 0180 
00190 
0 0200 
0 0210 
0 0220 
0 0230 
0 0240 
00250 
0 0260 
0 02  70 
0 0280 
0 0290 
0 0300 
0 0310 
0 0 320 
0 0330 
0 03  40 
00350 
0 0360 
0 0370 
0 0 380 
0 0390 
0 0 400 
0 0410 
0 0 420 
0 0430 
0 0 440 
0 0 450 
0 0460 
0 0470 
0 0 480 
0 0 49  0 
0 0 500 
00510 
0 0520 
00530 
0 0 540 
00550 
0 0 560 
0 0 570 
0 0 580 
0 0590 


PROGRAM  LIPDATECINPUT,OUTPUTj  TAPEl,  TAPE2>  TAPE3^ 

TAPE4^  TAPES,  TAPC:6,  TAPF  7,  TAPE8,  TAPE9  ) 

2 CONTINUE 

PRINT,  THIS  IS  DATA  PILE  UPDATE* 

PRINT,*  ENTER  FILE  NO* 

READ, IFILE 

PRINT,*ENTER,K  AND  KOUNT,  K=1  OR  2 FOR  FILE  TAPEl  OR  TAPE2  + 
>^RINT,  *KOUNT=NUMBER  OF  RECORDS  TO  BE  INSERTED* 

READ,K,KOUMT 

CALL  RETRC 5HTAPE1 ) 

CALL  RETHC2, 5HTAPE2) 

GO  TO  ( 1,  I,  1, 4,  5,  6,  7),  lElLE 

1 PRINT,  *ERROR- I FILE  LESS  THAN  4,  TRY  AGAIN* 

PRINT,  I FILE 
GO  TO  2 
4 NAMEs5HTA‘^E4 
GO  TO  88 
h NAME=5HTAPE5 
GO  TO  88 

6 NAME--5HTAPE6 
GO  TO  8 6 

7 NAME=5HTAPE7 
GO  TO  88 

8 NAME=5H TAPES 
GO  TO  88 

9 NAME-5HTAPE9 
GO  TO  88 

88  CALL  RETRC  1 FILE, NAME) 

REWIND  IFILE 

REWIND  1 

REWIND  2 

REWIND  3 

190  CONTINUE 

REAnciFiLF, 101)  ILAST 

PRINT, ILAST 

I EC ILAST. EG. 0)  GO  TO  888 
101  FORMATC I 10) 

NREC  = KOUNT+ILAST 
WRITEC3,101)  NREC 
PRINT, NREC 
I CO'JN=  ! 

130  CONTINUE 

READC I FILE,  1 00)  I, TIME 

100  FORMATC I 10, E20. 4) 

IFC ICOUN. EG. ILAST)  GO  TO  110 
WRITEC3,100)  I, TIME 
I COUN=  I COUN+  1 
GO  TO  130 

110  WRITEC3,100)  ICOUN, TIME 
I COUN=  ILAST+  1 
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0 0600 
00610 
00620 
00630 
0 0640 
00650 
00660 
00670 
00680 
0 0690 
0 0700 
0 0 710 
0 0 720 
0 0730 
0 0 740 
0 0750 
0 0760 
00770 
0 0780 
0 0 79  0 
0 0800 
00810 
00820 
00830 
0 0840 
00850 
00860 
00870 
00880 
00890 
00900 
00910 
0 0920 
0 0930 
0 09  40 
00950 
00960 
0 09  70 
0 0980 
0 0990 
0 1000 
0 1010 
0 1 020 
0 1030 
0 1 0 40 
0 1050 
0 1060 


150  KEWIND  K 
READCK, 101 ) I 

151  CONTINUE 
READCK,100)  I, TIME 
WRITE(3, 100)  ICOUN, TIME 

IF(  ICOUN.  EQ.NKFC)  GO  TO  140 

ICOUN=ICOUN-M 

GO  TO  151 

140  WRITEC3>100)  ICOUN^TIME 

FNDFILE  3 

REWIND  3 

REWIND  lEILE 

REAPC  3,  101 ) NREC 

WRITECIFILEj 101)  NREC 

240  READC3,  100)  1C0UN,TIME 

I FCICOUN. EO.NREC)  GO  TO  230 

WRI  TEC  I FILEj  100)  ICOIJN,  TIME 

GO  TO  240 

230  WRI  TEC  I FILE,  100)  I CO  UN,  TIME 

ENDFILE  IFILE 

REWIND  IFILE 

CALL  REPLC IFILE, NAME) 

PRINT, *FILE  UPDATE  IS, LAST  REC  CREATED  I 

PRINT,  IFILE,  ICOL«N 

GO  TO  99 

888  CONTINUE 

REWIND  IFILE 

NREC=KOUNT 

WRITECIFILE,  101)  NREC 
RFADCK,  725)  KOL>NT 
725  FORMATCI10) 

ICOUN=  1 

210  CONTINUE 

READC  K,  100)  I . TIME 

I FCI  .fQ.KOUNT)  GO  TO  200 

IFCICOUN.EQ.KOUNT)  GO  TO  200 

WRI  TEC  I FILE,  ’ 00)  ICOUN,TIME 

ICOUN=ICOL»')+l 

GO  TO  2 1 0 

200  WRI  TEC  IFILE,  100)  ICOUN,TlME 

PfjT  T . * r j L p UPDATE  ISsLAST  REC  CREATED'^' 

PRINT,  1 FILE,  I COUN 

FNDFILE  IFILE 

REWIND  IFILE 

CALL  REPLC IFILE, NAME) 

99  CONTINUE 
END 
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00100  PROGRAM  FILECRC INPUT>OUTPUT^ TAPF1» TAPF2, 

0 0110+  TAPE6, TAPE7^  TAPE8>  TAPE9 , TAPF 1 0. TAPFl 1,  TAPFlg) 
00100  DO  100  I = A.  10 

0 0130  GO  TO  C 4^  A,  4,  4,  5>  6.  7,  8,  9,  10»  1 1>  12),  I 

00140  4 NAME--5HTAPE1 

00150  IFII,E=1 

0 0160  GO  TO  110 

0 0170  5 \'AMF=5HTAPE0 

00180  1FILF=0 

0 0190  GO  TO  110 

00000  6 NAMF=5HTAPE6 

00210  1FIL£=6 

0 0220  GO  TO  !10 

00230  7 NAME=5HTAPF7 

00240  IFILF=7 

0 02  50  '0  TO  110 

00260  8 NAMF=5HTAPF8 

00270  IFILE-8 

0 0280  GO  TO  110 

00290  9 NAME=5HTAPE9 

00300  IFILFr!9 

0 0310  GO  TO  110 

00320  10  NAME=6HTAPF10 

00330  IF1LF=10 

0 0340  GO  TO  110 

0 0350  11  NAME--=6HTAPE11 

00360  IFILF=11 

0 0370  GO  TO  110 

00380  12  NAME=6HTAPE12 

00390  1F1LE=12 

00400  1 10  CONTINUE 

00410  REVIND  IFILE 

00420  ICOUN=0 

00430  11=1 

00440  TIMF=0. 

0 0450  WRI  TE<  I FILE,  30)  I COLIN 
00460  30  FORMAT!  I 10) 

00470  U'KI  TEC  IFILF,  35)  II, TIME 
0 0480  35  FORMAT!  I 10,  E20.  4) 

0 0490  FNDFILF  IFILE 
00500  REWIND  IFILE 
00510  100  CONTINUE 

0 0 520  FND 
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PROGRAM  LISTING  FOR  CONTINUOUS  SEARCH  - 
DECOY  FIEtJ)  SUBMODEL  (MOD  3) 


C-L 


0 03 


0 0 1 00 
0 0 110 
0 0 120 
00130 
0 0 140 

0 0 1 r>o 
00160 
0 01  70 
00180 
00190 
0 0200 
0 02  10 
0 0220 
0 0230 
0 0240 
0 02S0 
0 0260 
0 0 2 70 
0 020  0 
0 0 2 9 0 
0 0 30  0 
00310 
0 0320 
0 0330 
0 0 340 
0 0380 
0 0 3 60 
0 0 3 70 
0 0 38  0 
0 0 39  0 
0 0 400 
0 0 410 
0 0420 
0 0430 

0 0 4 40 

0 0 480 
0 0 460 
0 0 4 70 

O 0 4H0 

0 0 49  0 
0 0800 
0 0 81  0 
0 0 820 
0 0 830 
0 0540 
0 0 8 80 
0 0 560 
0 0 5 70 
0 0 5K0 
0 0 590 
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P^ROGRAM  M0D3C  I NPUT.  OUTHl'T.  lAPF  1.  TAPF2) 
n I MEN  b I ON  TMD(  5)  > XfK  81  . Yn<  ANOn(  8)  . INDC  5)  . CO  OKI)  ( 5) 

DIMENSION  NOFD<  5)  , I CHASIK  S)  , VnC(  51 
DIMENSION  TDCVA(  500)  . TOGO  Y ( 500  K PPH'  500) 

DIMENSION  VDC5) 

DIMENSION  DTM0N(5) 

DIMENSION  IPAR(5) 

PRINT,  ♦enter  NO  OE  DECOYS* 

RFAD,NDECOY 


DO  2000  NI=l,NDFCOY 
PRINT  2010, NI 

2010  E0RMAT<1H  , center  for  DFCOY  N0  + , I 3,  2X, +COURD,  VD*  ) 
HEAP, COURDCNI ) , VD(N I ) 

2000  CONTINUE 

PRINT,  +ENTER  VC,  VS,  COUKC,  COURS* 

HEAD,  VC,  VS,  CnukC.  COUkS 


PRINT,  +ENTEH  1C,  TM,  K 1 K,  K'LU-:* 

HEAD,  TC,  TM, K1 H, H2R 

PP^NT, +ENTER  NO  OF  REPL  I C/\TI  ON  SC  NO  T 10  EXCEED  500)* 

HEAn,NREPLI 

NINT=5O0 

I TGT5=0 

I PM  = 1 

IP2=5 


K=200. 

PRINT, *C* 

PRINT, *000,  DECOYS  KANDQm  START  IN  AREA4- 
PRINT, ♦SUB  RANDOM  START  ON  HOUNDARY* 
RRINT,*C* 

PRINT,  ♦COUNTER  TIME  FNCOl’NTERED+ 

CALL  SECONT'^X) 

Y=RANF(X) 

XX-0. 


TIMF=0. 


HRS=0. 


KO  IJNT=0 

’ ON 0=0 

K0NC=O 

15  CONTINUE 

no  662  !nY=1,NDFC0Y 

DTMONC I DY)=0, 

662  CONTINUE' 

T I M E=  0 . 

HRS=0. 

I N C=  0 
I N S=  0 
TNC=0. 

TN  S- 0 . 

I CHASS=  0 

DFr-=  87.  + C 1 7 . / 60.  ) 
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00600  no  0 I=l,NnF.COY 
0 0 6 10  ! N IX  I ) = 0 

0 0620  TNPC I ) = 0. 

00630  ICHASDCn  = 0 
0 0 640  2 CONTINUF 

00650  200  CONTINUF 

0 0 6A0  ANOC=  6.  283  19*KANF(  XX) 

0 0670  RADUS=h!*RANF<  XX  ) 

00680  XC=RADUS>i«COS(  ANPC) 

0G690  YC=RADIJS*SIN(  ANCC) 

0 0 700  TNC=TNC+1. 

00710  ANRS=  6.  28319+KANF(  XX) 

0 0 720  RAD=K 
0 0 730  X3=KAD+COSCANG83 
00740  YS=RAD*SIN(:ANRS) 

0 0 750  TNS=TNS+1. 

00760  INC=1 
0 0 770  INS=1 

00780  I FCNDECOY.  Ft).  0)  00  TP  800 

00790  DO  4 I=1,NPEC0Y 

00800  CALI.  I NnC0Y<  ANOIX  I ),  XLX  I ),  YD(  I )^  XC^  YC>  XS,  YS, 
0 06  10+  XX,  R,  VPC<  I )/ V'5,  VC) 

00820  TNPC  1 ) = TND(  I )+ 1 . 

0 08  30  IND(I)=l 
00840  4 CONTINUF 

00R50  800  CONTINUF 

0 0860  DT-0. 5 

00870  202  CONTINUE 

0 0880  TIMF- TI  MF+DT 
00890  H.RS=HR8+DT 
0 0900  IC  = 0 

0 0910  CALL  LITRAL<TNC,XC,YC,ANOC,XNOl,  YNOW,  INC, 

0 0920+  VC,  COURC,  XX,  VU,  V.S,  VDC,  X(),  YP,  I C,NPFCOY, 

0 0930+XC, YC) 

0 09  40  XC  = XNOW 
0 09  50  YC  = YNOW 

2 0960  1 F(  I CHA.5S.  FQ.  0)  00  TO  300 

P0970  1 F(  STMON.NF'Xfl'VlF)  00  TO  305 

0 0980  DO  310  I = 1,NDFC0Y 

0 0990  I FC  I CHASnC  r 3 . FO-  1 ) 00  TO  315 

01000  310  CONTINUF 

0 1010  PRINT,  + EKH~  FNCOUNTF.RF  D DFCOY  NOT  TURNED  0N  + 

0 1020  PRINT  3^ I 

3 1030  3 FORMATC  110) 

0 10  40  00  TO  00202 

0 1050  315  XS^XIKI) 

0 1060  Y5  = YD( I ) 

0 10  70  IN  5^1. 

01080  AN05=6.2R319*RANFCXX) 

0 1090  IN  5=1 
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01100  ICHASS»0 
01110  GO  TO  00320 
0 112fl  300  CONTINUE 

01130  IC*88 

0 1140  CALL  L1TBAL<TNS,  XS,YS,ANGS,XNO^',YNOV.v  IN.S,  VS/ 

0 1 150+  COURS,XX,VC/ VS>  VDCjXD/YD,  IC/NDECOY/ 

0 1 160+XC/YC) 

0 1170  XS-XNOW 
01180  YSsYNOW 
0 1190  GO  TO  00320 

0 1200  305  IF(TIME,LT.  STMONl  GO  TO  320 

0 1210  GO  TO  00300 

01220  320  CONTINUE 

0 1230  IFfNDECOY.EO.01  GO  TO  370 

0 1240  DO  330  I*l/NDECOY 

0 1250  IF(ICHASD(1).EQ.  n GO  TO  331 

01260  IF«TINE.LT.DTMON(I))  GO  TO  330 

0 1270  IF'TIME.EQ.DTMONn  1 ) GO  TO  330 

0 12B0  350  CONTINUE 

01290  IC=I 

0 1300  CALL  LITBAL<TND(I)/XD<I)/YD(n,ANGnU),XNOV'*YNOW> 
0 1310+  IND(  11/  VD(  I ),  COURIK  I 1 .XX,  VC,  VS,  VDC,  XD,  YD, 

0 1320+IC,NDECOY,XC,YO 
0 1330  XD(I)  = XNOW 
0 1340  YDUlsYNOW 
01350  GO  TO  00330 
0 1360  331  CONTINUE 

0 1370  I F<  TIMF.LT.  DTMON(  in  GO  TO  350 

0 1380  IFCTIME.EQ.DTMONdn  GO  TO  345 

01390  IF<( TIME-DTMON( I n .LT.DT)  GO  TO  342 
0 1400  PRINT,  ♦ERR- I CHASD=l,DTMON.  LT.  TIME* 

01410  PRINT,  TIME,  DTMONC  I ),  I 
0 1 420  I CHASDC I )=0 
01430  GO  TO  00330 
01440  342  CONTINUE 

0 1 450  I CHASDC I )=0 
01460  ICHASS=0 
p 1 /j70  DTMn\'<  n = DTMONC  n + TM 
0 1480  X.S=XDd) 

0 1490  YS=YDd) 

0 1500  IN  5=1 
0 1510  TNS=2. 

01520  ANGS=6.28319+RANE(XX) 

01530  GO  TO  00202 
01540  345  CONTINUE 

0 1550  ICHASDdn0 
0 1 560  DTMONd  ) = DTMONd  ) + TM 
01570  1CHASS=0 
01580  XS»XD(I) 

0 1 590  YSsYIH  I ■) 
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0 I 600  ANG?=6.08  319*1-ANF(X)<> 

01610  !NS= 1 
0 1600  TN.S-O. 

0 1630  r^O  TO  00  3 30 
0 1 F/iO  330  COF'TINHF 
0 16S0  3 70  CON  TIN OF 

0 1660  /iPO  CONTINUF 

0 1670  OALl..  MFFKVC.  VS.  XC,  YC,XS>  YF>  IO.NJ) 

0 1 68  0 IFdFN.FQ.n  00  TO  100 

0 1 6^50  00  TO  00390 

01700  390  CONTINUF 

01710  I F( I CKASS. FG. n 00  TO  00? 

0 1700  I FCNFIFCOY.  FC.  [,!)  00  TO  200 

0 173'^  ISAVF=0 

O17A0  no  /|/|0  I=l,NnFCOY 

0 17SO  I F(  1 CHAsn(  n . FG.  n oo  to  200 

0 1 760  I F<  TIMF.LF.  DTMONC  I ) ) OO  TO  AAVi 

01770  GALL  MFFT(  VnC<  n . VS.  XD(  I ) . Yl)(  I ) . XS,  Yb.  I IN) 

0 1 7B0  IF<1FN.F0.0)  00  TO  440 

0 1 790  nrSM=SGHT(  <xnc  I )-XS)*+.p+CYD(  I >- VS)t  + 2) 

01800  ISAVF=T 

01S1O  440  CONTINUF 

01820  I FC I SAVF. FO. 0)  00  TO  002 

0 1830  no  450  Isl.NL'FCOY 

OIF'140  I FC  I . FO'oI  SAVK)  00  TO  450 

01850  IF CTIMF.LF.DTMONC n ) oO  TO  450 

01860  GALL  MFFT C VUCf J ) . Vb. XD< I ). Y0( I ) , XS. Y S. 1 KN ) 

0 1H70  IF(IF,N.FG.0)  00  TO  450 

0 ! 8 80  ni  S=  SQRT(  ( Xl'H  I)  -XSY*  + 2+  < Y0(  I ) - YS)**0) 

01890  I F( ni S. OT. 01 SM)  00  TO  450 

01900  DISM=DIS 

0 19  10  ISAVF-I 

0 1920  450  CONTI, NUF 

0 1930  OTMONf  I ,SAVF)  = TIMF  *-TC 

0 1940  I CHASrU  I SAVF)  = 1 

0 1950  STMON=TIMF;-t  TC 

0 1 9 60  I CHAbb= 1 

01970  00  TO  00460 

0 10  8 0 100  KO  UN  T=  KO  UN  T+ 1 

0 1990  KOl\'C=kONC+ 1 

0 2000  TOf.VACKONCI -HKb 

0 2010  IPPINT^IPl';  1 

0 2020  I F(  I PK’INT.MF.  KOUNT)  00  TO  1000 

0 2030  IPl=f,  P1+' 

02040  PRINT  95.K0UNT.  rnCVACKONC) 

0PO50  95  FORMATf  1 1 0.  F 1 1 . 4> 

02060  KHITFC1.95)  KO  UN  T.  TDC  VA<  KONG  ) . T1  N F. 

0 2070  1000  CONTINUF 

0 2080  I FCKOUNT.  FU.MKFPLI  ) OO  TO  9 7 

02090  ISTAT-IP2+10 
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PIP100  IFfISTAT.NF.KOUNT)  C^O  TO  1010 
0P110  HALL  STATC  TDCVA>KONC>  FEANC>  STDFVC.  VAKC) 

0F1S0  PKINT>  *MEANC,  STDFVCj  VARf> 

0 P 1 30  PRINT, FRANC, STDFVC, VARC 

001^0  CALL  FHEOCTDCVA,KONC,  ISTAT.NRFPLI,NINT,  I TGTfi) 
0P150  PRINT, tFKFQUFNCY  GT  5000*-r  A+ 

0P160  PRINT, ITGT5 

0PI70  IFCKOND.EO.0)  GO  TO  15 

0P180  CALL  FRFQC  TDCOY,KONU,KOND,KOND,  NINT,  I TGT5) 
0P190  PRINT,  t-FRFGUFNCY  GT  5fl00.-DFCOY  + 

0P200  PRINT,  ITRT5 
0S210  IP2=IP2+5 

O2PP0  GO  TO  15 
02230  4^0  CONTINUF 

02240  KOND=KONn+l 
02250  TDC0YCK0ND)=HR5 
022R0  KPD(KOND)  = KOUNT+ 1 
0 22  70  GO  TO  00202 
02280  1010  CONTINUF 

0 2290  GO  TO  15 
0 2300  97  CONTINUF 

02310  CALL  STATI T0CVA,K0NC, FRANC, STDFVC, VARC) 

02320  PRINT, *MEANC, STAFVC.  VAKC* 

0 2330  PRINTl  76,  FFANC,  STOF.VC,  VARC 
0 2340  176  FORMAT?  3F1  1 . 4) 

02350  CALL  FRFQC  TDCVA,KONC,  I STAT,NRRPLI  ,NINT,  I TGT5) 
02360  PRINT, ♦FRFQUFNCY  GT  5000.-CVA+ 

0 2370  PRINT,  ITGT5 
0 2380  GO  TO  999 

O23'’0  I FfKONn.  FO.  0)  GO  TO  999 

02400  CALL  FKFOCTDCOY,KOND,KONn,KOND.NINT,  I TGT5) 
02410  PRINT, *FRFQUFNCY  GT  5000. - DPCOY+ 

0 2420  PRINT,  ITGT5 

0 2430  I F<kONH.  RP.  0)  GO  TO  999 

G 2440  CALI,  .STAT<  TDCOY  . KOND,  FFANI),  .STOFVD,  VARl) ) 

0 2 45M  PR  IN 'I,  + Ml-  AND,  STDEVD,  VARD* 
p 2460  PRINT  1 76.,  FFAND,  .STDEVD,  VAK'D 
p 999  CONTINUF 

0 2480  1 I 1 CONTINUF 

0 2 490  FND 
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0 PSRP! 
0 P5  1 R 
n ?SPR 
R P53R 
R PS/jO 
R P5S0 
055(«,0 
F)  PS7fl 
0 P5BR 
0 P 59  R 
R P ^ H f) 

M p ,<  I 3 

R 26PR 
0 P^3R 
0 2 6/10 


niMFN.SI^Ow''  STDFV,  VAK) 

•S!!M=0, 

•^0  3 Isi.icoiiT 

-■■’.'Ms  SUM+  TI  MFf  n 
^ CONTINUE 
EEANr  SUN/r  LOATH  COUT) 

VAR=0. 

DO  10  j_  j 
yAR=VAKH(TlMF(n.pF.^^,,, 

«0  CONTINUE 

VAR=VAK/EL0ATnC0l'T> 

stofv=sgktcvako 

kp'Turn 

FND 


L, 


..Lj-iU 


LE  COPY 
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0 2 650  5UBK0UTINF  LI  TBALf  TN,  XLAST»  YLAST/  ALA5T,  XNOV,  YNO  W>  IND, 
02660+  \;,COUKSF#XX,VC,  VS,  VDC,Xn,YD,  IC,NL'FCOY, 

02670+  XC,YC> 

02680  DIMFN'SION  DYMRC  5)  , SUMRC  5) 

02690  DIMENSION  Xnf  , Y[U' 5)  , ' OCC  5) 

0 2700  DT=.5 

0 2710  K=200. 

0 2 720  IF(IC.EU.0)  RO  TO  20 
02730  IF(  in.rO.BB)  RO  TO  30 

0 27A0  I F(  I C.LE.  NDFCOY)  GO  TO  25 

02750  DK’INT,*FKR>  rC  INCORKECT+ 

0 2760  PRINT,  I C 

0 2770  RETURN 

0^780  20  CONTINUE 

0 2 790  IFC1NU.F&.0)  GO  TO  15 

0C3R0O  CALL  XYNOWC  TN,  PT,  V,  XLAST,YLAST,  COURSE,  ALAS1,  K, 

0 28  1 0+  XNOV;,  YNOW,  IND,  ANOV,  1 ANRI 
02820  65  DO  m I=l,NOE'COY 

02830  DYMR(I)  = .SQKTC  (XD(  I ) -XNO  '.■0  + + 2+  f YDC  I ) - YNO  W)  + + 2> 

0 28  60  RCS=R2R 
0 28  50  Vt'--VDC<I) 

02860  SlIMRin -RCS+Rl  K 

02870  IFCDYMPd  1 .LE.  SUMRf  I >)  RO  TO  35 

0 2«80  60  CONTINIIF 

0 2890  IFnANR.FO,0)  GO  TO  999 

02900  ALAST=ANOV 

0 2910  IANR=0 

0 2920  GO  TO  9 99 

0 22  30  3 5 IN!')=0 

0 2960  15  CALL  NRH  Y ( XX,  R,  V,  "T,  XL  AST,  YL  AS  l Al.Ab’l,  AND  W, 

0 29  50+  XNO  V,  YNO  , i AN  R) 

0 29  60  TN=2. 

0 29  70  lNn=i 
0 2980  IANG=1 
099(7,  AO  T O 6 5 
0 3000  25  CONTINUE 

0 3010  I E ( I N n . F (■: . 0 ) f-0  T 0 50 

0 3020.  CALL  XYNOVT  TN,  DT,  V,  XL.ASI  , YLAST,  COURSE,  ALAST,  R, 

0 3030+  XNOl-  , YNOW,  IND,  ANOV,',  I ANG) 

0 30^0  RC5-K2.P 

0 3 0 50  60  DYMRC I C)  = SUKT<  C XC- XNO V) + + 2+ ( . C- YNO  W) + + 2) 

0 3060  VR=VDC(IC) 

0 3070  SIIMRC  I C)  = r:CS+ I'll  R 

03080  I EC  PYMRC  I C)  . LF.  SUMK<  I C)  ) GO  TO  55 
0 3fJ9O  I EC  I ANR- EG.  0)  RO  TO  999 
0 3100  ALAST=AN0W 
0 3110  IANG-0 
0 3 120  GO  TO  9 99 
0 3130  55  COM  TIN!  IF 

0 3160  IND=0 
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0 31  SO 
0 3 160 
0 3 I 70  + 
0 3 1 B 0 
0 3190 
0 3200 
0 321  0 
0 3220 
0 3230 
0 32^0 
( 3250+ 
0 32  60 
0 32  70 
0 3280 
0 329  0 
0 3300 
0 3310 
0 3320+ 
0 3 3 30 
0 3340 
0 3350 
0 3360 
0 33  70 
0 3380 
n 3390 


50  CONTINUE 

IND=  1 
TN=2. 

I ANn= 1 
CO  TO  60 
30  CONTINUE 

IF<INn.EO.0)  RO  TO  65 

lEfIANG.EO.0)  po  TO  999 
ALAST=ANOW 
I ANG=0 
GO  TO  999 
65  CONTINUE 

XNOr  VNnu/J 
XNOv.S  rNOW#IANA> 

TN=2. 

IND-I 


ALASTsANOV 

I AN  P=  0 

999  CONTINUE 

RETURN 

END 
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0 3A0a 
0 34  J0 
0 34«0 
0 3430 
0 3440 
0 34S0 
0 3 4A0 
0 3 470 
0 3480 
0 3 49  0 


0A/1 1 /71 . 


SUPHOUTI  NF  MFF.K  VC*  VS»  .XC*  YC*  XS*  YS*  I EN  ) 
K0=60. 

RANGF-SGKrC ( X C~ X 8) * * ?+ ( YC- Y S) * ♦ 2 ) 

I FC  RAl\)GF.  FT.  H0)  CO  TO  10 
I FN=  1 

GO  TO  00020 
10  IFN=0 
20  CONTINUE 
RFTURN 
FND 
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0 3600 
0 3610  + 
0 3620 
0 3530 
0 3540 
0 3550 
0 3560 
0 3 570 
0 3580 
0 3590 
0 3600 
0 3610 
0 3 620 
0 3 630 
0 3640 
0 3650 
0 36  60 


SUBROUTINF  INDCOY(ANGDl..Xni,Yni^XC.YC,  XS,  YS,  XX,  K, 
VDC,  VS>  VO 
5 CONTI  NUF 
AN  on  1 = 6. 28  3 19+ RAN F( XX) 

RAD=R+RAN1-'(XX) 

XDl=RAn+C05(  ANGDl  ) 

Y01=RAD+SINC ANCDl ) 

CALL  MEETC  VDC,  vs,  XDl,  YDl,  XS,  YS,  I EN) 

IFaEN.FO.l)  on  TO  5 

RnS=Rl H 

RCS=R2R 

SUM  R=RDS+RCS 

RCVAnY=SQRT( ( XC- XD 1) * + 2+ ( Y(>  YD  1 ) * + 2) 

I RCVAPY.  FT,  SUMR)  RETURN 

GO  TO  5 

RETURN 

cND 
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0 SL'RRrMrriNF  XYMOWC  TN^  DT>  V>  XLA.ST>  YLA.^T,  COIJKiiF^  ALAST^ 

B 3A8B+  XNQ  V,  YNO  V,  IND,  ANO  R,  I ANA) 

B 369A  CONTI  NUF 

A3  70P.  FI  ,FN=TN*DT*  V 

0 371C  1 F<  FLFN.  AT.  COUKSF.)  GO  TO  10 

037?n  nX=cnT*V)+C0S(ALA5T) 

0 3739  DY- (OT+V)*SIN(  ALAST) 

0 37A0  RN0V.'=SQKT<(DX+XLAST)**2+(DY+YLAST)+  + 2) 

0 3750  IFCRNOW.GT.  R)  GO  TO  10 
R37A0  XNOW=XLAST+ DX 
0 3770  YN0V.'=YLAST+DY 
0 3 780  TN=TM+i. 

0 3790  IANG=0 
0 3800  GO  TO  99 
A3E10  10  ( ONTINUF 

0 3»20  CALL  NPHYf  XX,  R,  V,  DT,  XLAST,  YLAST,  ALAS7,  ANOlv,  XNO  F , YNOF 
0 38  30  TN=g. 

0 38/(0  INn=l 
0 3B50  IANG-1 
0 3KX0  99  CONTINUF 
0 3870  RFTIIPN 
0 9880  FNri 


1 ANG) 
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n 7,p,9P 
P 3900  + 
0 3910 

( ,i  39C50 

0 39  30 
0 39  4 0 
0 39  bO 
0 39  60 
0 3970 
0 3980 
0 3990 
0 4000 
0 40  10 
0 40P0 
0 40  30 


S'irkODTINF  NPHY(XX,  V^DT^XLAST.  YLAi;T,  ALAST, 
ANO  VP  XNO  VP  YNO  VP  1 AN  0) 

20  CONTINUE 
ALPHA=6.28319*KANF ( XX) 

PHY=ALPHA+ALAST 
PKY= AMOfK  PHY,  A. 283 1 9 ) 
nX=  ( DT*  V)  >t=COS(PHY) 

DY=C  DT*V) *SIN(PHY) 

RNOV.'=SQRT(  ( DX+XLAST)**2+  <DY  + YLAST>**2) 

I F( RNOV. GT. R)  GO  TO  20 

ANOV'=PHY 

XNOW=XLAST+DX 

YN0V'=YLA.9T+DY 

RFTUHN 

END 
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p/10/iO  SUBROL'TiMF  FHEQC  TUCVAj  KOi\'C»  I STAT,  NKFPL  I . Nl  I NT.  ITGT5) 
0^050  DIMENSION  TDCVA( S0O) . I FHE< 5P0I . TCOM C 500) 

0/1060  DIMENSION  IF1^EC(  500) 

0A0  70  ITRT5=0 

0 A08O  IFC  ISTAT.lt. NRFPLI  ) 00  TO  5 

0 A090  ISLIM  = NREPLI 
0 4 100  00  TO  6 

0 4110  5 ISUM=ISTAT 

0 4120  6 TINT=10. 

0 4130  TMAX=0. 

0 41  40  DO  7 J=1.NINT 
0 4150  lFRF(J)-0 
0 4160  IFRFC(J)  = 0 
0 4170  7 CONTINUE 

0 4180  DO  10  1=1,1  SUM 

0 4190  DO  20  J=  NINT 


0 4200 
0 42  1 0 

TCOMPA=FLO,-.  r?  JJ  + TINT 
IF(  TDCVA<  I ) . GT.  5000.  ) 

GO 

TO 

40 

0 4220 

IF<  TDCVA< I ) . GT, TCOMPA) 

GO 

TO 

20 

0 42  30 
0 42  40 

IFRF(J)=IFRF(J)+1 
IFCTCOMPA.LT.  TMAX)  GO 

TO 

10 

0 42  50 

TMAX=TC0MPA 

0 4260 

GO 

TO  10 

0 42  70 

40 

I TGT5=1 TGT5+ 1 

0 4280 

GO 

TO  10 

0 429  0 

20 

CONTINUE 

0 4300 

10 

CONTTN'UF; 

0 43  1 0 

I P = 

TMAX/TINT4  1 . 

0 43  20  DO  30  1=1,  IP 

0 4330  TOOM( I ) =FLOAT< I 1 + 1 1 NT 
0 43  40  30  CONTINUE 

0 4350  I FHFCC  1 ) = I FKF<  1 > 

0 4360  DO  110  1 = 2,  IP 

O 43  70  I FKFC ( r ) = I FRF( I)  + I FhFC< I - 1 ) 

0 4"*RO  110  CONTINUE 

0 4390  PRINT, +TIME  I N TE  RVAL , FPEOUFNC  E,  Cl  W F* 
0 4400  DO  100  1=1,  IP 

04410  1 FC I FKFC I ) . FU. 0)  GO  10  100 

0 4420  PRINT  3 5, TCOMC I ) , I FRF ( I) , I FRFC( I ) 

0 4430  3 5 FORMAT  ( IH  , 2X , E 11  • 2,  2X,  21  1 0 ) 

0 4 440  100  CONTI N U E 

0 4450  RETURN 

0 4440  fND 
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